Generalizable links between borderline personality traits and functional connectivity.

BACKGROUND: Symptoms of borderline personality disorder (BPD) often manifest in adolescence, yet the underlying relationship between these debilitating symptoms and the development of functional brain networks is not well understood. Here we aimed to investigate how multivariate patterns of functional connectivity are associated with borderline personality traits in large samples of young adults and adolescents. METHODS: We used functional Magnetic Resonance Imaging (fMRI) data from young adults and adolescents from the Human Connectome Project: Young Adults (HCP-YA; N=870, ages 22-37 years, 457 female) and Development (HCP-D; N=223, ages 16-21 years, 121 female). A previously validated BPD proxy score was derived from the NEO Five Factor Inventory (NEO-FFI). A ridge regression model with cross-validation and nested hyperparameter tuning was trained and tested in HCP-YA to predict BPD scores in unseen data from regional functional connectivity. The trained model was further tested on data from HCP-D without further tuning. Finally, we tested how the connectivity patterns associated with BPD aligned with age-related changes in connectivity. RESULTS: Multivariate functional connectivity patterns significantly predicted out-of-sample BPD scores in unseen data in young adults (HCP-YA; pperm=0.001) and older adolescents (HCP-D; pperm=0.001). Regional predictive capacity was heterogeneous; the most predictive regions were found in functional systems relevant for emotion regulation and executive function, including the ventral attention network. Finally, regional functional connectivity patterns that predicted BPD cores aligned with those associated with development in youth. CONCLUSION: Individual differences in functional connectivity in developmentally-sensitive regions are associated with borderline personality traits.

The brain's topographical organization shapes dynamic interaction patterns that support flexible behaviour based on rules and long term knowledge.

Adaptive behavior relies both on specific rules that vary across situations and stable long-term knowledge gained from experience. The frontoparietal control network (FPCN) is implicated in the brain's ability to balance these different influences on action. Here, we investigate how the topographical organization of the cortex supports behavioral flexibility within the FPCN. Functional properties of this network might reflect its juxtaposition between the dorsal attention network (DAN) and the default mode network (DMN), two large-scale systems implicated in top-down attention and memory-guided cognition, respectively. Our study tests whether subnetworks of FPCN are topographically proximal to the DAN and the DMN, respectively, and how these topographical differences relate to functional differences: the proximity of each subnetwork is anticipated to play a pivotal role in generating distinct cognitive modes relevant to working memory and long-term memory. We show that FPCN subsystems share multiple anatomical and functional similarities with their neighboring systems (DAN and DMN) and that this topographic architecture supports distinct interaction patterns that give rise to different patterns of functional behavior. The FPCN acts as a unified system when long-term knowledge supports behavior but becomes segregated into discrete subsystems with different patterns of interaction when long term memory is less relevant. In this way, our study suggests that the topographic organization of the FPCN, as well as the connections it forms with distant regions of cortex, are important influences on how this system supports flexible behavior.Significance Statement Adaptive behavior depends on adjudicating between specific rules that vary across situations. The frontoparietal control network (FPCN) helps guide this process through its interactions with other brain regions. We examined how local topographical features support this function of the FPCN. Subnetworks within the FPCN share key anatomical and functional features with adjacent systems linked to external attention and long-term knowledge. This topographic architecture supports the emergence of distinct interaction patterns: FPCN subnetworks act cohesively when long-term memory can support behavior, but segregate when long-term memory is not aligned with current goals. Our study shows that, in addition to dynamic interaction with spatially distant cortical regions, local topographical features of the FPCN play a significant role in flexible behavior.

Diffusion MRI head motion correction methods are highly accurate but impacted by denoising and sampling scheme.

Head motion correction is particularly challenging in diffusion-weighted MRI (dMRI) scans due to the dramatic changes in image contrast at different gradient strengths and directions. Head motion correction is typically performed using a Gaussian Process model implemented in FSL's Eddy. Recently, the 3dSHORE-based SHORELine method was introduced that does not require shell-based acquisitions, but it has not been previously benchmarked. Here we perform a comprehensive evaluation of both methods on realistic simulations of a software fiber phantom that provides known ground-truth head motion. We demonstrate that both methods perform remarkably well, but that performance can be impacted by sampling scheme and the extent of head motion and the denoising strategy applied before head motion correction. Furthermore, we find Eddy benefits from denoising the data first with MP-PCA. In sum, we provide the most extensive known benchmarking of dMRI head motion correction, together with extensive simulation data and a reproducible workflow. PRACTITIONER POINTS: Both Eddy and SHORELine head motion correction methods performed quite well on a large variety of simulated data. Denoising with MP-PCA can improve head motion correction performance when Eddy is used. SHORELine effectively corrects motion in non-shelled diffusion spectrum imaging data.

XCP-D: A Robust Pipeline for the post-processing of fMRI data.

Functional neuroimaging is an essential tool for neuroscience research. Pre-processing pipelines produce standardized, minimally pre-processed data to support a range of potential analyses. However, post-processing is not similarly standardized. While several options for post-processing exist, they tend not to support output from disparate pre-processing pipelines, may have limited documentation, and may not follow BIDS best practices. Here we present XCP-D, which presents a solution to these issues. XCP-D is a collaborative effort between PennLINC at the University of Pennsylvania and the DCAN lab at the University at Minnesota. XCP-D uses an open development model on GitHub and incorporates continuous integration testing; it is distributed as a Docker container or Singularity image. XCP-D generates denoised BOLD images and functional derivatives from resting-state data in either NifTI or CIFTI files, following pre-processing with fMRIPrep, HCP, and ABCD-BIDS pipelines. Even prior to its official release, XCP-D has been downloaded >3,000 times from DockerHub. Together, XCP-D facilitates robust, scalable, and reproducible post-processing of fMRI data.

Network connectivity and local transcriptomic vulnerability underpin cortical atrophy progression in Parkinson's disease.

Parkinson's disease pathology is hypothesized to spread through the brain via axonal connections between regions and is further modulated by local vulnerabilities within those regions. The resulting changes to brain morphology have previously been demonstrated in both prodromal and de novo Parkinson's disease patients. However, it remains unclear whether the pattern of atrophy progression in Parkinson's disease over time is similarly explained by network-based spreading and local vulnerability. We address this gap by mapping the trajectory of cortical atrophy rates in a large, multi-centre cohort of Parkinson's disease patients and relate this atrophy progression pattern to network architecture and gene expression profiles. Across 4-year follow-up visits, increased atrophy rates were observed in posterior, temporal, and superior frontal cortices. We demonstrated that this progression pattern was shaped by network connectivity. Regional atrophy rates were strongly related to atrophy rates across structurally and functionally connected regions. We also found that atrophy progression was associated with specific gene expression profiles. The genes whose spatial distribution in the brain was most related to atrophy rate were those enriched for mitochondrial and metabolic function. Taken together, our findings demonstrate that both global and local brain features influence vulnerability to neurodegeneration in Parkinson's disease.

Functional connectome fingerprinting across the lifespan.

Systematic changes have been observed in the functional architecture of the human brain with advancing age. However, functional connectivity (FC) is also a powerful feature to detect unique "connectome fingerprints," allowing identification of individuals among their peers. Although fingerprinting has been robustly observed in samples of young adults, the reliability of this approach has not been demonstrated across the lifespan. We applied the fingerprinting framework to the Cambridge Centre for Ageing and Neuroscience cohort (n = 483 aged 18 to 89 years). We found that individuals are "fingerprintable" (i.e., identifiable) across independent functional MRI scans throughout the lifespan. We observed a U-shape distribution in the strength of "self-identifiability" (within-individual correlation across modalities), and "others-identifiability" (between-individual correlation across modalities), with a decrease from early adulthood into middle age, before improving in older age. FC edges contributing to self-identifiability were not restricted to specific brain networks and were different between individuals across the lifespan sample. Self-identifiability was additionally associated with regional brain volume. These findings indicate that individual participant-level identification is preserved across the lifespan despite the fact that its components are changing nonlinearly.

Neurophysiological signatures of cortical micro-architecture.

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6800 time-series features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is co-localized with multiple micro-architectural features, including gene expression gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Generalizable links between symptoms of borderline personality disorder and functional connectivity.

Background |: Symptoms of borderline personality disorder (BPD) often manifest in adolescence, yet the underlying relationship between these debilitating symptoms and the development of functional brain networks is not well understood. Here we aimed to investigate how multivariate patterns of functional connectivity are associated with symptoms of BPD in a large sample of young adults and adolescents. Methods |: We used high-quality functional Magnetic Resonance Imaging (fMRI) data from young adults from the Human Connectome Project: Young Adults (HCP-YA; N = 870, ages 22-37 years, 457 female) and youth from the Human Connectome Project: Development (HCP-D; N = 223, age range 16-21 years, 121 female). A previously validated BPD proxy score was derived from the NEO Five Factor Inventory (NEO-FFI). A ridge regression model with 10-fold cross-validation and nested hyperparameter tuning was trained and tested in HCP-YA to predict BPD scores in unseen data from regional functional connectivity, while controlling for in-scanner motion, age, and sex. The trained model was further tested on data from HCP-D without further tuning. Finally, we tested how the connectivity patterns associated with BPD aligned with age-related changes in connectivity. Results |: Multivariate functional connectivity patterns significantly predicted out-of-sample BPD proxy scores in unseen data in both young adults (HCP-YA; pperm = 0.001) and older adolescents (HCP-D; pperm = 0.001). Predictive capacity of regions was heterogeneous; the most predictive regions were found in functional systems relevant for emotion regulation and executive function, including the ventral attention network. Finally, regional functional connectivity patterns that predicted BPD proxy scores aligned with those associated with development in youth. Conclusion |: Individual differences in functional connectivity in developmentally-sensitive regions are associated with the symptoms of BPD.

Integrating multimodal and multiscale connectivity blueprints of the human cerebral cortex in health and disease.

The brain is composed of disparate neural populations that communicate and interact with one another. Although fiber bundles, similarities in molecular architecture, and synchronized neural activity all reflect how brain regions potentially interact with one another, a comprehensive study of how all these interregional relationships jointly reflect brain structure and function remains missing. Here, we systematically integrate 7 multimodal, multiscale types of interregional similarity ("connectivity modes") derived from gene expression, neurotransmitter receptor density, cellular morphology, glucose metabolism, haemodynamic activity, and electrophysiology in humans. We first show that for all connectivity modes, feature similarity decreases with distance and increases when regions are structurally connected. Next, we show that connectivity modes exhibit unique and diverse connection patterns, hub profiles, spatial gradients, and modular organization. Throughout, we observe a consistent primacy of molecular connectivity modes-namely correlated gene expression and receptor similarity-that map onto multiple phenomena, including the rich club and patterns of abnormal cortical thickness across 13 neurological, psychiatric, and neurodevelopmental disorders. Finally, to construct a single multimodal wiring map of the human cortex, we fuse all 7 connectivity modes and show that the fused network maps onto major organizational features of the cortex including structural connectivity, intrinsic functional networks, and cytoarchitectonic classes. Altogether, this work contributes to the integrative study of interregional relationships in the human cerebral cortex.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks.

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Individual differences in delay discounting are associated with dorsal prefrontal cortex connectivity in children, adolescents, and adults.

Delay discounting is a measure of impulsive choice relevant in adolescence as it predicts many real-life outcomes, including obesity and academic achievement. However, resting-state functional networks underlying individual differences in delay discounting during youth remain incompletely described. Here we investigate the association between multivariate patterns of functional connectivity and individual differences in impulsive choice in a large sample of children, adolescents, and adults. A total of 293 participants (9-23 years) completed a delay discounting task and underwent 3T resting-state fMRI. A connectome-wide analysis using multivariate distance-based matrix regression was used to examine whole-brain relationships between delay discounting and functional connectivity. These analyses revealed that individual differences in delay discounting were associated with patterns of connectivity emanating from the left dorsal prefrontal cortex, a default mode network hub. Greater delay discounting was associated with greater functional connectivity between the dorsal prefrontal cortex and other default mode network regions, but reduced connectivity with regions in the dorsal and ventral attention networks. These results suggest delay discounting in children, adolescents, and adults is associated with individual differences in relationships both within the default mode network and between the default mode and networks involved in attentional and cognitive control.

Development of top-down cortical propagations in youth.

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

Neurophysiological signatures of cortical micro-architecture.

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6 800 timeseries features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is colocalized with multiple micro-architectural features, including genomic gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Network structure and transcriptomic vulnerability shape atrophy in frontotemporal dementia.

Connections among brain regions allow pathological perturbations to spread from a single source region to multiple regions. Patterns of neurodegeneration in multiple diseases, including behavioural variant of frontotemporal dementia (bvFTD), resemble the large-scale functional systems, but how bvFTD-related atrophy patterns relate to structural network organization remains unknown. Here we investigate whether neurodegeneration patterns in sporadic and genetic bvFTD are conditioned by connectome architecture. Regional atrophy patterns were estimated in both genetic bvFTD (75 patients, 247 controls) and sporadic bvFTD (70 patients, 123 controls). First, we identified distributed atrophy patterns in bvFTD, mainly targeting areas associated with the limbic intrinsic network and insular cytoarchitectonic class. Regional atrophy was significantly correlated with atrophy of structurally- and functionally-connected neighbours, demonstrating that network structure shapes atrophy patterns. The anterior insula was identified as the predominant group epicentre of brain atrophy using data-driven and simulation-based methods, with some secondary regions in frontal ventromedial and antero-medial temporal areas. We found that FTD-related genes, namely C9orf72 and TARDBP, confer local transcriptomic vulnerability to the disease, modulating the propagation of pathology through the connectome. Collectively, our results demonstrate that atrophy patterns in sporadic and genetic bvFTD are jointly shaped by global connectome architecture and local transcriptomic vulnerability, providing an explanation as to how heterogenous pathological entities can lead to the same clinical syndrome.

Mapping neurotransmitter systems to the structural and functional organization of the human neocortex.

Neurotransmitter receptors support the propagation of signals in the human brain. How receptor systems are situated within macro-scale neuroanatomy and how they shape emergent function remain poorly understood, and there exists no comprehensive atlas of receptors. Here we collate positron emission tomography data from more than 1,200 healthy individuals to construct a whole-brain three-dimensional normative atlas of 19 receptors and transporters across nine different neurotransmitter systems. We found that receptor profiles align with structural connectivity and mediate function, including neurophysiological oscillatory dynamics and resting-state hemodynamic functional connectivity. Using the Neurosynth cognitive atlas, we uncovered a topographic gradient of overlapping receptor distributions that separates extrinsic and intrinsic psychological processes. Finally, we found both expected and novel associations between receptor distributions and cortical abnormality patterns across 13 disorders. We replicated all findings in an independently collected autoradiography dataset. This work demonstrates how chemoarchitecture shapes brain structure and function, providing a new direction for studying multi-scale brain organization.

neuromaps: structural and functional interpretation of brain maps.

Imaging technologies are increasingly used to generate high-resolution reference maps of brain structure and function. Comparing experimentally generated maps to these reference maps facilitates cross-disciplinary scientific discovery. Although recent data sharing initiatives increase the accessibility of brain maps, data are often shared in disparate coordinate systems, precluding systematic and accurate comparisons. Here we introduce neuromaps, a toolbox for accessing, transforming and analyzing structural and functional brain annotations. We implement functionalities for generating high-quality transformations between four standard coordinate systems. The toolbox includes curated reference maps and biological ontologies of the human brain, such as molecular, microstructural, electrophysiological, developmental and functional ontologies. Robust quantitative assessment of map-to-map similarity is enabled via a suite of spatial autocorrelation-preserving null models. neuromaps combines open-access data with transparent functionality for standardizing and comparing brain maps, providing a systematic workflow for comprehensive structural and functional annotation enrichment analysis of the human brain.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex.

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Local molecular and global connectomic contributions to cross-disorder cortical abnormalities.

Numerous brain disorders demonstrate structural brain abnormalities, which are thought to arise from molecular perturbations or connectome miswiring. The unique and shared contributions of these molecular and connectomic vulnerabilities to brain disorders remain unknown, and has yet to be studied in a single multi-disorder framework. Using MRI morphometry from the ENIGMA consortium, we construct maps of cortical abnormalities for thirteen neurodevelopmental, neurological, and psychiatric disorders from N = 21,000 participants and N = 26,000 controls, collected using a harmonised processing protocol. We systematically compare cortical maps to multiple micro-architectural measures, including gene expression, neurotransmitter density, metabolism, and myelination (molecular vulnerability), as well as global connectomic measures including number of connections, centrality, and connection diversity (connectomic vulnerability). We find a relationship between molecular vulnerability and white-matter architecture that drives cortical disorder profiles. Local attributes, particularly neurotransmitter receptor profiles, constitute the best predictors of both disorder-specific cortical morphology and cross-disorder similarity. Finally, we find that cross-disorder abnormalities are consistently subtended by a small subset of network epicentres in bilateral sensory-motor, inferior temporal lobe, precuneus, and superior parietal cortex. Collectively, our results highlight how local molecular attributes and global connectivity jointly shape cross-disorder cortical abnormalities.

Brain atrophy progression in Parkinson's disease is shaped by connectivity and local vulnerability.

Brain atrophy has been reported in the early stages of Parkinson's disease, but there have been few longitudinal studies. How intrinsic properties of the brain, such as anatomical connectivity, local cell-type distribution and gene expression combine to determine the pattern of disease progression also remains unknown. One hypothesis proposes that the disease stems from prion-like propagation of misfolded alpha-synuclein via the connectome that might cause varying degrees of tissue damage based on local properties. Here, we used MRI data from the Parkinson Progression Markers Initiative to map the progression of brain atrophy over 1, 2 and 4 years compared with baseline. We derived atrophy maps for four time points using deformation-based morphometry applied to T1-weighted MRI from 120 de novo Parkinson's disease patients, 74 of whom had imaging at all four time points (50 Men: 24 Women) and 157 healthy control participants (115 Men: 42 Women). In order to determine factors that may influence neurodegeneration, we related atrophy progression to brain structural and functional connectivity, cell-type expression and gene ontology enrichment analyses. After regressing out the expected age and sex effects associated with normal ageing, we found that atrophy significantly progressed over 2 and 4 years in the caudate, nucleus accumbens, hippocampus and posterior cortical regions. This progression was shaped by both structural and functional brain connectivity. Also, the progression of atrophy was more pronounced in regions with a higher expression of genes related to synapses and was inversely related to the prevalence of oligodendrocytes and endothelial cells. In sum, we demonstrate that the progression of atrophy in Parkinson's disease is in line with the prion-like propagation hypothesis of alpha-synuclein and provide evidence that synapses may be especially vulnerable to synucleinopathy. In addition to identifying vulnerable brain regions, this study reveals different factors that may be implicated in the neurotoxic mechanisms leading to progression in Parkinson's disease. All brain maps generated here are available on request.

Multimodal phenotypic axes of Parkinson's disease.

Individuals with Parkinson's disease present with a complex clinical phenotype, encompassing sleep, motor, cognitive, and affective disturbances. However, characterizations of PD are typically made for the "average" patient, ignoring patient heterogeneity and obscuring important individual differences. Modern large-scale data sharing efforts provide a unique opportunity to precisely investigate individual patient characteristics, but there exists no analytic framework for comprehensively integrating data modalities. Here we apply an unsupervised learning method-similarity network fusion-to objectively integrate MRI morphometry, dopamine active transporter binding, protein assays, and clinical measurements from n = 186 individuals with de novo Parkinson's disease from the Parkinson's Progression Markers Initiative. We show that multimodal fusion captures inter-dependencies among data modalities that would otherwise be overlooked by field standard techniques like data concatenation. We then examine how patient subgroups derived from the fused data map onto clinical phenotypes, and how neuroimaging data is critical to this delineation. Finally, we identify a compact set of phenotypic axes that span the patient population, demonstrating that this continuous, low-dimensional projection of individual patients presents a more parsimonious representation of heterogeneity in the sample compared to discrete biotypes. Altogether, these findings showcase the potential of similarity network fusion for combining multimodal data in heterogeneous patient populations.